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Abstract 



Parametric simultaneous solitary wave (simulton) excitations are shown pos- 
sible in nonlinear lattices. Taking a one-dimensional diatomic lattice with 
a cubic potential as an example we consider the nonlinear coupling between 
the upper cutoff mode of acoustic branch (as a fundamental wave) and the 
upper cutoff mode of optical branch (as a second harmonic wave) . Based 
on a quasi-discreteness approach the Karamzin-Sukhorukov equations for two 
slowly varying amplitudes of the fundamental and the second harmonic waves 
in the lattice are derived when the condition of second harmonic generation 
is satisfied. The lattice simulton solutions are explicitly given and the results 
show that these lattice simultons can be nonpropagating when the wave vec- 
tors of the fundamental wave and the second harmonic waves are exactly at 
ix I a (where a is the lattice constant) and zero, respectively. 

PACS numbers: 63.20. Pw, 63.20. Ry 
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I. INTRODUCTION 



Since the pioneering work of Fermi, Pasta, and Ulam[l] on the nonlinear dynamics in lat- 
tices, the understanding of the dynamical localization in ordered, spatially extended discrete 
systems have experienced considerable progress. In particular, the lattice solitons, which 
are localized nonlinear excitations due to the balance between nonlinearity and dispersion 
of the system, are shown to exist, and many important applications are found in transport 
of energy, proton contactivity, structural phase transition and associated central-peak phe- 
nomena, etc [2,3]. In recent years, the interest in localized excitations in nonlinear lattices 
has been renewed due to the identification of a new type of anharmonic localized modes[4-6]. 
These modes, called the intrinsic localized modes (ILM's)[4], or discrete breathers [5, 6], are 
the discrete analog of the lattice envelope (or breather) solitons with their spatial extension 
only a few lattice spacing and the vibrating frequency above the upper cutoff of phonon 
spectrum band. The ILM's have been observed in a number of experiments [7- 13]. Recently, 
much attention has been paid to the gap solitons in diatomic lattices[14-25]. In linear case, 
a diatomic lattice allows two phonon bands. There is a upper cutoff for phonon frequency 
and a frequency gap (forbidden band) between acoustic and optical bands, induced by mass 
and/or force-constant difference of two different types of particles. No interaction occurs 
between phonons and the phonons can not propagate in the system when their frequencies 
are in the gap or above the phonon bands. However, these properties of the phonons change 
drastically when nonlinearity is introduced into the system. New types of localized modes, 
especially the gap solitons, may appear as nonlinear localized excitations with their vibrat- 
ing frequencies in the band gap. The gap solitons and ILM's as well as their chaotic motion 
have been observed in damped and parametrically excited one-dimensional (ID) diatomic 
pendulum lattices [26-28]. 

On the other hand, in recent years numerous achievements have been made for optical 
solitons in nonlinear optical media[29-31]. Besides the temporal optical solitons, which is 
promising for long-distance information transmission in fiber, spatial optical solitons also 
attract much attention. The spatial optical solitons are believed to be the candidates for 
all-optical devices, such as optical switches and logic gates, etc[32]. Recently, the study 
of optical parametric processes, in particular the second harmonic generation (SHG), which 
marked the birth of nonlinear optics, has generated a great deal of new interest [33]. It 
was suggested that it is possible to obtain large nonlinear phase shifts by using a cascaded 
second-order nonlinearity [34]. In 1974, Karamzin and Sukhorukov (KS) recognized that 
the cascaded second-order parametric processes may support solitons under general phase- 
matching conditions. They derived two coupled nonlinear equations for the envelopes of 
the fundamental and second harmonic waves [35]. The difference between the KS equations 
and the envelope equations for usual SHG is the inclusion of dispersion and/or diffraction, 
which are necessary for short pulses and/or narrow light beams. Simultaneous solitons (i. e. 
two components are solitons) are found for the KS equations and these solitons are later 
termed as the parametric simultons[36]. The concept of the simultons has been generalized 
to the nonlinear optical media with periodically varying refractive index [3 7]. Since the 
eigenspectrum of linear electromagnetic waves consists of many photonic bands and the 
vibrating frequencies of the simultons may be in the gaps between these bands, the name 
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parametric band-gap simulton is given by Drummond et al[37-39]. Different from the self- 
trapping mechanism for Keer solitons, the formation of the simultons is due to the energy 
transfer and mutual self-trapping between the fundamental and the second harmonic waves. 

On the contrary, the SHG in lattices is much less investigated. Although in the standard 
textbook of solid state physics[40] there exists a simple experimental description for three 
phonon processes in solids, it seems that there is no detailed theoretical approach to the 
SHG in nonlinear lattices until recently. In a recent paper, Konotop considered theoretically 
the SHG in a nonlinear diatomic lattice and obtained some interesting results[41]. 

In many aspects a nonlinear lattice is similar to a nonlinear, periodic optical media. The 
discreteness of lattice results in the symmetry breaking for continuous translation and makes 
the property of the system periodic, in particular the frequency spectrum of corresponding 
linear wave splits into many bands. It should be stressed that the SHG does not occur 
in ID monatomic lattices (see the next section). However, a SHG can be realized if we 
consider nonlinear multi-atomic lattices . The reason is that in the monatomic lattices, an 
efficient energy transfer (resonance) between any two modes in the system does not occur. 
But the situation is different for the multi-atomic lattices. A multi-atomic lattice allows 
many branches of linear dispersion relation, and the dispersion relation is periodic with 
respect to lattice wave vector. It is just the multiple- value and periodic property of the 
dispersion relation makes it possible that the phase-matching condition for the SHG, i.e. 
the condition by which the resonance between the fundamental and second harmonic waves 
may occur, can be satisfied by selecting the wave vectors and the corresponding frequencies 
from different spectrum branches. 

Motivated by the study of the optical simultons, in this paper we show that lattice 
simultons are possible in the multiatomic lattices with cubic nonlinearity (different from the 
case in nonlinear optics, here the order of nonlinearity means the order in the Hamiltonian 
of the system). The paper is organized as follows. The next section presents our model and 
an asymptotic expansion based on a quasi-discreteness approach. In section III we solve the 
KS equations derived in section II and provide some lattice simulton solutions. A discussion 
and summary is given in the last section. 



II. MODEL AND ASYMPTOTIC EXPANSION 



A. The model 



As mentioned in the last section, the SHG may in princile occur in any multi-atomic 
lattice, but for definiteness and for the sake of simplicity we consider here a ID diatomic 
lattice with a cubic interaction potential. The Hamiltonian of the system is given by 




+-K 3 (w n - v n f + -K' 3 (v n+1 - w n ) 3 + ±V 3 vl + -Vlwl 



(1) 
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where v n —v n (t) (w n = w n (t) ) is the displacement from its equilibrium position of the nth 
particle with mass m (M). n is the index of the nth unit cell with a lattice constant a = 2a®, 
ao is the equilibrium lattice spacing between two adjacent particles. Here for generality we 
assume that the nearest-neighbor force constants Kj(j = 2, 3) in the same cells are different 
from the nearest-neighbor force constants Kj(j = 2,3) in different cells. V 3 and V 3 ' are the 
force constants related to the on-site cubic potential for two types of particles. Without loss 
of generality we assume m < M, Kj < Kj(j = 2, 3), and V 3 ' < Vj$. The equations of motion 
for describing the lattice read 

d 2 

-^V n = I 2 (w n - V n ) + l' 2 {w n -\ ~ V n ) + h(w n - V n f - l' 3 (w n -i ~ V n f - Q! m U^, (2) 

d 2 

— w n = J 2 (v n - w n ) + J 2 (v n+1 - w n ) - J 3 (v n - w n ) 2 + J' z (v n+l - w n ) 2 - a M wl, (3) 

where Ij = Kj/m,^ = K'-jm,Jj = KjjM,J'- = K'-jM (j = 2,3), a m = V 3 /m and a M = 
VI I 'M. The linear dispersion relation of Eqs. (2) and (3) is given by 

u±(q) = ^{(h + I 2 + J2 + J 2 ) ± [{h + I' 2 + J 2 + J' 2 f ~ !QhJ 2 sin 2 (ga/2)] "} 2 , (4) 

where the minus (plus) sign corresponds to acoustic (optical) mode. Thus we have two 
phonon bands u±(q) and obviously u±(q + Q) = oo±(q), here Q = 2jir/a, j is an integer and 
Q is the reciprocal lattice vector of the system. At the wave number q — 0, the phonon 
spectrum has a lower cutoff u;_(0) = for the acoustic mode and an upper cutoff c^+(0) = 
(I 2 + I 2 + J 2 + J' 2 ) l l 2 for the optical mode. At q = ir/a there exists a frequency gap between 
the upper cutoff of the acoustic branch u^(n/a) and the lower cutoff of the optical branch 
u + (Ti/a), where w±(n/a) = (l/v^){(/ 2 + V 2 + J 2 + J' 2 ) ± [(/ 2 + V 2 + J 2 + J' 2 f - 16/ 2 ^] 1/2 } 1/2 . 
The width of the frequency gap is uj + {ti /a) —u^{ji/a), which approaches zero when m — > M 
and K' 2 — * K 2 . This is just the limit of monatomic lattice with the lattice constant ao = a/2. 
We assume the gap is not small, i. e. we have (1 — m/M) and (1 — K 2 /K' 2 ) are of order unity. 

Because of the periodic property of u±(q), the condition of a second harmonic resonance 
in the system (2) and (3) reads 

g 2 = 2 Ql + Q, (5) 
uj 2 = 2u u (6) 



where q\ (q 2 ) and u\ (cu 2 ) are the wave vector and frequency of the corresponding fundamen- 
tal (second harmonic) wave, respectively. Eqs. (5) and (6) are also called the phase-matching 
conditions for the SHG. It is easy to show that in the limit m — > M and K' 2 — > K 2 the 
condition (5) and (6) can not be satisfied except for zero-frequency mode, i. e. the SHG is 
impossible in monatomic lattices. For the diatomic lattice, in order to fulfil (5) and (6) we 
may chose cu 1 G 0J-{q) and u 2 G u> + (q), then the conditions (5) and (6) give 

\l 2 + I' 2 + J 2 + J' 2 ) 2 - AI 2 J' 2 sin( gi a)l 1/2 
= 3(/ 2 + I' 2 + J 2 + J' 2 ) - 4 



(I 2 + I' 2 + J 2 + J' 2 ) 2 - m 2 J' 2 sin 2 (g!a/2) 



1/2 



(7) 
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It is available to solve q\ from the above equation. For simplicity we consider cutoff modes 
of the system. We take q\ = n/a,q 2 = and Q = —2ir/a, then the condition (5) is 
automatically satisfied. The condition (6) (the same as (7) ) now reads 

I 2 + I' 2 + J 2 + J' 2 = ^=^hA. (8) 

Eq.(8) also means that uj x = cj_(vr/a) = (1/2)(J 2 + V 2 + J 2 + J' 2 ) 1 ' 2 = (AI 2 J' 2 /3) 1 ^ and 
uj 2 = uj + (0) = (7 2 + I' 2 + J 2 + J' 2 y/ 2 = 2(4/ 2 J^/3) 1/4 If the all harmonic force constants are 
equal, i.e. K' 2 = K 2 , Eq.(8) gives m = M/3. Another particular case is all masses are the 
same, i.e. m = M. In this case Eq.(8) requires K 2 = K 2 /3. In general, the phase-matching 
conditions (5) and (6) impose a constraint on masses and harmonic force constants of the 
lattice. 



B. Asymptotic expansion 



We employ the quasi-discreteness approach (QDA) developed in Refs. 17 and 24 for 
diatomic lattices to investigate the SHG in the system (2) and (3). We are interested in the 
cascading processes of the system in which the width of excitation is narrower than usual 
SHG case. Thus we use different spatial-temporal scales in deriving the envelope equations 
for the fundamental and the second harmonic waves. We make the expansion 

^(t) = e [ U W + e l ' 2 u^ n + eu (2 l + •••], (9) 

where u n {t) represents v n (t) or w n (t). e is a smallness and ordering parameter denoting the 
relative amplitude of the excitation and u^ n = u^(^ n ,r; ra (t)), with 

£ n = e 1 ' 2 (na-\t), (10) 
r = et, (11) 
<p n = qna - u(q)t, (12) 

where A is a parameter to be determined by a solvability condition (see below) . Substituting 
(9)-(12) into Eqs.(2) and (3) and equating the coefficients of the same powers of e, we obtain 
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- l2(w® - v<£l) - J^gU - vjj} n ) = M n % (13) 
M n °l = 0, (14) 
M« - 2A ^ v m - T'n—w {0) M ^ 

M W - o\ ^ - ( 2— 4 A 2 — ^ «(°) 4- V (-n—w {1) 4- — —w {0) ) 

Mn > n ~ 2X dtdC n ' n [ 2 dtdr +X dej n ' n 2 1 dt™^' 1 + 2! dH n ' n ~ 1 ) 

+h(w$> n ~ <i) 2 - />i°Li - v^f - a m (vZ) 2 , (16) 
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- Uvgn ' ~ *Sh-i - w®) = N®, (17) 



and 

N%1 = 0, (18) 

N ^ = 2X ^k~n W ^ + ^Wn^ (19) 

AT(2) - 2A ^ W W - ( 2— + X 2 —) w<® + J' (a—v {1) , + —— v (0) ,\ 

-J"3«i - ^°i) 2 + - «ffl 2 - «"«l) 2 , (20) 



with j = 0, 1, 2, • • •. Eqs.(13) and (17) can be rewritten in the following form 

Lw^l = J i M nl + + (J^ + h + I^j Ni% (21) 

[dt 2 

where the operator L is defined by 

= (j^ + h + I'^j (jj + -h + J^j u^ n - {I 2 J 2 + I' 2 J' 2 



+ h + tfi = /awgj, + + M{1 (22) 



(?) 
n,n 

-/ 2 J£ + «2U), (23) 

where u^ n (j = 0, 1, 2, • • •) are a set of arbitrary functions. From Eq.(21) we can solve w$ n 
and obtain a series of solvability conditions (envelope equations) whereas Eq.(22) is used to 

C. Envelope equations for cascading processes 

We now solve Eqs.(22) and (23) order by order. For j — it is easy to get 



= Ax(r, Cn) exp(^) + A 2 (r, £ n ) exp(*0+) + c. c, (24) 
S = -T-TTr Al ( T ' ex P(^« ) + '2 T \ T » ^(r, a) exp(^0+) + c. c. (25) 



with 0^ = gna — w± (<?)£. cf±(g) have been given in Eq.(4). A 1 and A 2 are yet to be 
determined two envelope (or amplitude) functions of the acoustic and the optical excitations, 
respectively. They are the functions of the slow variables £ n and r. c. c. denotes the 
corresponding complex conjugate. For simplicity we specify two modes, i. e. the acoustic 
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upper cutoff mode (q 1 = n/a, cu 1 = uj^{ji/a) = {AI^J^/?)) 1 ^ ) and the optical upper cutoff 
mode (q 2 = 0, to 2 = u+(0) = 2uo x = 2(4 J 2 </ 2 /3) 1/4 )• Thus we have 

= Ai(r, £n)(-l) n exp(-iwit) + A 2 (r, f n ) exp(-iw 2 f) + c. c, (26) 
«SS = , >2 7i ' " \ r, Ai(r,en)(-l) n exp(-M0 

— + i 2 + i 2 
I + I' 

+ 2 2 , T 2 , t , A 2 (t, f n ) exp(-^ 2 t) + c. c. (27) 

— w 2 + i 2 + i 2 

From the discussion in subsection II. A, the modes chosed in such way satisfy the phase- 
matching conditions (5) and (6) for the SHG. Thus in Eqs.(27) and (28) A\ (A 2 ) represents 
the amplitude of the fundamental (second harmonic) wave, respectively. 

In the next order (j=l), a solvability condition of Eqs.(21) and (22) requires the pa- 
rameter A = 0, thus £ n = na. The second-order solution reads 

w£l = B + [Bi(-l) n exp(-Mt) + B 2 exp(-iu 2 t) + c. c] , (28) 
<!n = B +\ (-1) exp(-zwit) 

(h + I' 2 )B 2 - r 2 adA 2 /d(j n . \ 

+ ,2,r —J, exp(-^ 2 t) + c. c. , (29) 

— L0 2 + 1 2 + 1 2 J 

where Bj (j=0, 1,2) are undetermined functions of £ n and r. 

In the order j—2, we have the third-order approximate equatuion 

LwW = J 2 M$ + J' 2 M^ +l +(J^ + I 2 + I'^j Ni% (30) 

Eq.(22) is not necessary since from (30) we can obtain closed equations for A 1 and A 2 . Using 
Eqs.(26)-(29) we can get M^l, M^ n+1 and N^ n . By a detailed calculation we obtain the 
sovability condition of Eq.(30) 

: dAi 1 (PAi 

dr 2 1 9^ 

.QAx 1 d 2 A 2 

dr + 2 2 9a 



^ + -r 1 ^ + A 1 ^A 2 = 0, (31) 
i^ + ^r 2 ^ + A 2 A 2 = 0, (32) 



where the coefficients are expressed as 

2 = "^[-A,- 1 - A 2 (/ 2 + J 2 )(J 2 + ' ( } 

a [1 - A 2 (/ 2 + I 2 )]A 3 - ArW - AiA 2 (/ 2 2 - (/^) 2 )(J 2 - J> m , . 

A = A 4 - \ 2 l a M - A 2 (I 2 - J 2 ) 2 (J 2 + J 2 )a M (3Q) 

2 2cu 2 [A 2 r 1 + A 2 (/ 2 + / 2 )(J 2 + J^)] ' 1 J 
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with 

U 2 \ ■ „ , 0' = 1 ) 2) (37) 

—a; j + i 2 + J 2 

A 3 = (/s - /£)( J 2 - 4) - K\Js - J' 3 ) + (h - I' 2 ) [(J 3 + JS) - Ai(/ 3 + J£)(J 2 - J' 2 )\ , (38) 
A 4 = [1 - Ax(/ 2 - / 2 )] 2 [-J 3 A 2 1 + / S (J 2 + 4)\ + [1 + Ai(/ 2 - mV's^ 1 ~ W2 + (39) 



Introducing the transformation Uj = eAj (j = 1, 2) and noting that £ n = e 1 / 2 x n (i n = no) 
and r = ei, Eqs.(31) and (32) can be rewritten into the form 

* + ^ + Al „,^ , (40, 

^ + ^ + A2llN 0. (4!) 



We should point out that Eqs.(5) and (6) are perfect phase-matching conditions for the 
SHG. If we allow a small mismatch for frequency, Su, the conditions (5) and (6) become 

cu 2 = 2wi + 6u, q 2 = 2gi + Q. (42) 



In this case Eqs.(40) and (41) change into 



' dui dui\ 1 d 2 ui . „ . . 

."a +2r 1 ^ r + A ltI ; tl2 exp(- lfot )=0, (43) 

du 2 du 2 \ 1^ <9 2 w 2 .9 /.c. x , ,, s 



where i>j (j = 1, 2) are the group velocities of of the fundamental and the second harmonic 
waves near at q = ir/a and q = 0, respectively. 

Eqs.(43) and (44) are the coupled-mode equations for the fundamental and the second 
harmonic waves. Such equations have been obtained by Karamzin and Sukhorukov in the 
context of nonlinear optics [35]. One of important features of the KS equations is the inclusion 
of dispersion, which is absent in usual SHG envelope equations [41]. 



III. LATTICE SIMULTON SOLUTIONS 



In this section, we solve the KS equations (43) and (44) derived in our lattice model 
and thus present some lattice simulton solutions for the system (2) and (3). In general, the 
property of the solutions of Eqs.(43) and (44) depends strongly on the coefficients appearing 
in the equations, in particular on their signs. At first we notice that in our system, I\ and 
T 2 , which are respectively the group- velocity dispersion of the fundamental and the second 
harmonic waves, are both negative. But the signs of the nonlinear coefficients, A 1 and A 2 , 



8 



may be generally of both signs. Thus the situation here is different from the KS equations 
derived for the cascading process in nonlinear optics, where the nonlinear coefficients have 
the same sign, while the group- velocity dispersions may have different signs [42]. 

To solve Eqs.(43) and (44), we make the transformation 



ui{x n) t) = E/i(C)exp[i(A;ia; n - Slit)], 
u 2 (x n , t) = U 2 (C) exp[i(k 2 x n - fl 2 t)}, 



(45) 
(46) 



with ( = kx n — Qt. Substituting (45) and (46) into (43) and (44), we obtain 

d 2 XJ x 



d( 2 
d 2 U 2 

d( 2 



+ a x U x U 2 - fi x U x = 0, 



+ a 2 U 2 - p 2 U 2 = 0, 



(47) 
(48) 



where a x = 2A l /(T l k 2 ),a 2 = 2A 2 /(T 2 k 2 ), f3 l = -2(Q X - Vl k - \T .k 2 ) / (T 1 k 2 ) , (3 2 = -2(Q 2 - 
v 2 k — ^T 2 kl)/ (T 2 k 2 ), Q = v\k + Y\kk\, with k 2 = 2ki,Q 2 = 2Qi + 5uj and k\ = (v 2 — 
vi)/(Ti — 2T 2 ). One of the coupled soliton-soliton (i. e. simultaneous solitons for two wave 
components) solutions of Eqs. (47) and (48) reads 



^OL\OL 2 VO 



where a condition (5\ — (5 2 — —4 is required. The parameter k is given by 

2{y 2 - Vl )h + (fi - 2T 2 )k 2 + 5uj 



(49) 
(50) 



k 



2(r 2 - 2i\) 



(51) 



The lattice configuration in this case takes the form 



w n (t) = (-iy 

12 

CCi 

Vn(t) = (-1)" 

12 



12 



— sech 2 (kna — Vtt) 



\JOL\a 2 
2 

- — sech 2 (kna — Qt) 



cos[kina — {uj\ + £li)t] 



12 



y/®ia 2 -uj( + I 2 + I 2 



cos[k 2 na — (uj 2 + fl 2 )t], 
-2 



sech (kna — £lt) 



cos[kina — (uo\ + £l\)t\ 



h + r* 



«i -u\ + i 2 + r 2 



- — sech 2 (kna — Qt) 



cos[k 2 na — (u 2 + Q 2 )t}. 



(52) 



(53) 



If qi (q 2 ) is exactly equal to 7i/a (zero) but with Su ^ 0, one has v± — v 2 — 0. In this case 
k x = k 2 = 0, fix = 2Y x k 2 ,Sl 2 = 2T 2 k 2 ,n = and k = {5lu/[2(T 2 - 2F!)]} 1 / 2 . (52) and 
(53) present a nonpropagating simulton excitation, in which the vibrating frequency of the 
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acoustic- (optical-) mode component being within the acoustic (optical) phonon band. In our 
model, T 2 — 2Fi > thus 5oo should be taken positive in this case. In addition, from (52) 
and (53) we see that the envelopes for both the acoustic and optical components are kinks 
(or dark solitons). Futhermore, if K' 2 = K 2 , the displacement of light particles, v n (t), only 
has an optical-mode component. 

The other simulton solution of Eqs.(47) and (48) reads 

C/i = ■ sechV, (54) 

6/4 



U 2 = - sech^C , (55) 

CKi \0 



where we have f3i = —f3 2 = —4. The parameter k now reads 



_ 2{v 2 - vi)ki + 2(r 2 - T 1 )k 2 1 - 5u; 

k ~ 2(2^ + r 2 ) • (56) 



The lattice configuration is now given by 



12 

w n {t) = (— l) n+1 sech 2 (kna — tit) cos[k\na — (ui + tii)t] 

yJOL\OL 2 

( sech 2 (A;na — tit)) cosf^na — {uj 2 + ti 2 )t], 

«i V3 / 



(57) 



v n (t) = (-1) 



n+l 



12 



h-I' 2 



■ k — sech (kna — tit) cosf/cma — (u\ + tii)t\ 

Jaiai-(4 + h + r 2 



12 i 2 + r 2 

ai -u\ + i 2 + r 2 



sech 2 (kna — tit) 



cos[k 2 na — (u 2 + ti 2 )t}. 



(58) 



Thus in this case the acoustic-mode component is a staggered envelope soliton but the 
optical-mode component is still an envelope kink. If v± — v 2 — we have k± — k 2 — 
0,ti 1 = 2T 1 k 2 ,ti 2 = -2T 2 k 2 ,ti = and k = {-5lu/[2(T 2 - 2F1)]} 1 / 2 . In this situation the 
simulton (57) and (58) is also a nonpropagating excitation with the vibrating frequency of 
the acoustic- (optical-) mode component within(above) the acoustic (optical) phonon band. 
In order to make k to be real we should take 5uo < in this case. 

A common requirement for the existence of the simulton solutions (52), (53), (57) and 
(58) is sgn^a^) > 0, which means sgn(A!A 2 ) > because > in our model. It can 
be met by chosing different values of system parameters. For example, in the following two 
particular cases we have sgn(AiA 2 ) > 0: 



1. K' 2 = K 2 , K' 3 = K 3 = 0. In this case Ai = -a M /ui, A 2 = - J 2 a M /[2^ 2 (/ 2 + J 2 )}. 

2. K' 2 = K 2 , VI = V 3 = 0. In this case Ax = (J' 3 - J 3 )(l + I 2 /J 2 )/uj u A 2 = (I' 3 - J 3 + 
4-Js)/[2u; 2 (I 2 + J 2 )]. 
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IV. DISCUSSION AND SUMMARY 



We have analytically shown that the lattice simultons are possible in nonlinear diatomic 
lattices. Based on the QDA for the nonlinear excitations in diatomic lattices developed 
before [17, 24], we have considered the resonant coupling between two phonon modes, one 
from the acoustic and other one from the optical branch, respectively. The KS equations 
are derived for the envelopes of the fundamental and second harmonic waves by taking 
new multiple spatial-temporal scale variables, which are necessary for narrower nonlinear 
excitations. Exact coupled soliton (simulton) solutions are obtained for the KS equations 
and the simulton configurations for the lattice displacements are explicitly given. 

Similar to the optical simultons in nonlinear optical media, the physical mechanism 
for the formation of the lattice simultons is due to the cascading effect between two lattice 
wave components. In this process, the fundamental and the second harmonic waves interact 
with themselves through repeated wave-wave interactions. For instance the energy of the 
fundamental wave is first upconverted to the second harmonic wave and then downconverted 
again, resulting in a mutual self-trapping of each wave thus the formation of two simultaneous 
solitons. 

Mathematically, in addition to the resonance conditions (5) and (6), the formation of 
a lattice simulton needs a balance between the cubic nonlinearity (in the Hamiltonian) and 
the dispersion, the latter is provided by the discreteness of the system. Thus for deriving the 
envelope equations in this case, we must chose the multiple-scale variables different from the 
ones used in usual SHG. In our derivation for the KS equations based on the QDA [17,24], 
only one small parameter, i. e. the amplitude of the excitation, is used. This method gives a 
clear, justified and self-consistent hierarchy of scales and thus the corresponding solvability 
conditions, which are just the envelope equations we need. Thus it is satisfactory according 
to the point of view of singular perturbation theory. 

Cubic nonlinearity exists in most of realistic atomic potentials [24]. Thus it is possible 
to observe the lattice simultons reported here. It must be emphasized that the multi-value 
property of the linear dispersion relation is important for generating the simultons in lattices. 
Thus a diatomic or multi-atomic lattice is necessary for observing such excitations. 

The theory given above can be applied to multi-atomic and higher-dimensional lat- 
tices, and higher-order nonlinearity can also be included. For instance, if we consider the 
Hamiltonian with cubic and quartic nonlinearities, Eqs.(31) and (32) sould be generalized 
to 

r)A 1 r) 2 A 

i-j^- + + Ai A\A 2 + (AnlAI 2 + A l2 \A 2 \ 2 )A 1 = 0, (59) 

*TF + \ V ^W + A2 A ' + (A2l|A|2 + A ^ A ^ A2 = °> ( 6 °) 

where Aij(i,j = 1,2) are self-phase and cross-phase modulational coefficients contributed 
by the quartic nonlinearity of the system. Eqs.(59) and (60) can be derived using the 
multiple-scale variables £ n = ex n ,r = e 2 t under the assumption v n (t) = 0(e), w n (t) = 0(e), 
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K 3 = 0(K' 3 ) = 0(e), and V 3 = 0(V£) = 0(e). A small frequency mismatch can also be 
included in (59) and (60) and similar equations like (43) and (44) with additional self- and 
cross-phase modulational terms can also be written down. A detailed study will be presented 
elsewhere. 
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